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Abstract 

In this paper we will review a recently introduced method for solving the 
Hamilton-Jacobi equations by the method of Separation of Variables. This 
method is based on the notion of pencil of Poisson brackets and on the bi- 
hamiltonian approach to integrable systems. We will discuss how separability 
conditions can be intrinsically characterized within such a geometrical set-up, 
the definition of the separation coordinates being encompassed in the bihamil- 
tonian structure itself. We finally discuss these constructions studying in details 
a particular example, based on a generalization of the classical Toda Lattice. 
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1 Introduction 



The study of the separability of the Hamilton-Jacobi (H J) equations associated with a 
given Hamiltonian function if is a very classical issue in Mechanics, dating back to the 
foundational works of Jacobi, Stackel, Levi-Civita and others. It has recently received 
a strong renewed interests thanks to its applications to the theory of integrable PDEs 
of KdV type (namely, the theory of finite gap integration) and to the theory of 
quantum integrable systems (see, e.g., [TT] IHlj). 

As it is well known, the problem can be formulated as follows. Let (M, u) be a 2n 
dimensional symplectic manifold, and let (pi, . . . ,p n , q%, ■ ■ ■ , q n ) = (p, q) be canonical 
coordinates (in U C M), i.e., U\ v = Ym=i ^Vi A dq{. 

The (stationary) HJ equation reads 

H{q 1 ,...,q n ;— ,...,—) = E . (1.1) 

OPl OPn 



Definition 1.1 A complete integral 5(q; a%, . . . , a n ) of HJ is a solution of 
depending on n parameters (a±, . . . , a n ) such that Det d ® 2 g a ^ 
H is said to be separable in the coordinates (p, q) if HJ admits an additively separated 
complete integral, that is, a complete integral of the form 

n 

5(q; a 1 ,...,a n ) = ) j S^; ax,..., a n ). (1.2) 

i=i 

In this paper we will focus on an equivalent 'Definition-Theorem' of separability, 
originally due to Jacobi and recently widely used by Sklyanin and his collaborators. 
Let us consider an integrable Hamiltonian H, that is, let us suppose that, along with 
H = H\ we have further n — 1 mutually commuting integrals of the motion H2, . . . H n , 
with dHi A ... A dH n ^ 0. 

Definition 1.2 An integrable system (Hi,...,H n ) is separable in the coordinates 
(p, q) if there exist n non-trivial relations 

$ i (q i ,p i ;H 1 ,...,H n ) = Q, i = l,...n, (1.3) 

connecting single pairs {qi,Pi) of canonical coordinates with the n Hamiltonians Hi. 
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We called this a Definition-Theorem due to the fact that it is a constructive approach 
to separability, since the knowledge of the separation relations (|1.3|) allows one to 
reduce the problem of finding a separated solution of HJ to quadratures. Indeed, 
one can solve the relations $j(gj,pj; Hi, . . . , H n ) = 0, i — 1, . . . n for the pi to get 
Pi — Pi(qi', Hi, . . . , H n ) and then define: 

S(q\a u ...,a n ) = 22 Pi(li, H x , . . . , H n ) i dq[ . (1.4) 

i=i J ' H i= a i 

This is by construction a separated solution of HJ; the fact that it is a complete 
integral is equivalent to the fact that the integrals of the motion depend non trivially 
on the momenta. 

In intrinsic terms, one notices that, the equations Hi = ati, % = 1, . . . , n define 
a foliation T of M. The leaves of T are nothing but the (generalized) tori of the 
Arnol'd-Liouville theorem. The foliation is Lagrangian, that is, the restriction of 
the two-form uj to T vanishes. Hence the restriction to T of the Liouville form 
9 = YH=iPid<li to T is (locally) exact. Indeed, the function S defined by (jl.4j) is 
a (local) potential for the restriction of 6 to T . What is non intrinsic, and singles 
out the separation coordinates (p, q) is that the separation relations ()Oj) . which are 
another set of defining equations for the foliation T , have the very special property of 
containing a single pair of canonical coordinates at a time. The problem to find such 
a system of coordinates and relations is the core of the theory of SoV. In particular, 
a natural question arises: 

Is it possible to formulate intrinsic condition(s) on the Hamiltonians (Hi, . . . , H n ) 
to a priori ensure separability in a set of canonical coordinates? 

Actually, this is the main issue studied by both the 'classical' Eisenhart-Benenti 
theory of separability of natural systems defined on cotangent bundles to Riemannian 
manifolds (M,g), as well as the 'modern' theory, due to the St. Petersburg and 
Montreal schools, of SoV for systems admitting a Lax representation. We notice 
that both such general approaches require the presence of an additional structure 
to solve the problem. Indeed, the Eisenhart-Benenti theory requires the existence 
of a conformal Killing tensor for the metric g, while the Lax theory requires - in 
addition to the knowledge of a Lax representation with spectral parameter for the 
Hamiltonian system under study - the existence of an r-matrix structure for such a 
Lax representation. 

The method we will present /review in this paper, which has recently been exposed 
in the literature (see, e.g., [Zl E01 E H31 El 1201 [^1 EH E2I UHl IU ) can be seen as a kind 
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of bridge between the classical and the modern points of view, putting an emphasis 
on the geometrical structures of the Hamiltonian theory. Its 'additional' structure 
is simply the requirement of the existence, on the symplectic manifold (M, to) of a 
second Hamiltonian structure, compatible with the one defined by u. Namely, the 
bihamiltonian structure on M will allow us: 

1. To encompass the definition of a special set of coordinates, to be called Darboux- 
Nijenhuis (DN) coordinates, within a well defined geometrical object. 

2. To formulate intrinsic (i.e., tensorial) conditions for the separability of a Hamil- 
tonian integrable system, in the DN coordinates associated with the bihamil- 
tonian structure. 

3. To give recipes to characterize, find and handle sets of DN coordinates. 

The paper is organized as follows. In Section |21 we will briefly introduce the notions 
of bihamiltonian geometry relevant for the paper. In particular, we will discuss the 
notion of DN coordinates, as well as methods to find them In Section El we will 
discuss the main theorems of the bihamiltonian set-up for SoV, namely, the tensorial 
conditions ensuring Separability of the HJ equations in DN coordinates. 

Then we will discuss our constructions in a specific example, whose separability, 
to the best of our knowledge, has not yet been considered in the literature. It is a 
generalization of the periodic Toda lattice with four sites. In Section H] we will recall 
its definition, and show how the "bihamiltonian recipe" for SoV can be applied to it. 
Although our constructions can be generalized to the generic iV-site system, for the 
sake of concreteness and brevity we choose to consider the four-site system only, and 
sometimes rely on direct computations to prove some of its properties. 

In the last subsection we will apply our geometrical scheme to study a specific 
reduction of this generalized Toda system, and find integrals of the motion which are 
not encompassed in the Lax representation. This result can be possibly as first step 
towards an alternative approach to the so-called chopping method of [S] for the full 
(non-periodic) Toda Lattice. 
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2 Some issues in the geometry of (bi)hamiltonian 
manifolds 



We start this Section recalling some well known facts in the theory of Poisson mani- 
folds. 

Definition 2.1 A Poisson manifold (M, {•, •}) is a manifold endowed with a Poisson 
bracket, that is a bilinear antisymmetric composition laws defined on the space C°°(M) 
satisfying: 

1. The Leibnitz rule: {fg, h} = f{g, h} + g{f, h}; 

2. The Jacobi identity {/, {g, h}} + {g, {h, /}} + {h, {/, g}} = 0. 

We recall that a Poisson bracket (or Poisson structure) can be equivalently described 
with a Poisson tensor, that is, with an application P : T*M — > TM, smoothly varying 
with m G M defined by: 

{f,g} = (df,Pdg) 

where (•, •) denotes the canonical pairing between T*M and TM. In a given coor- 
dinate system (x 1 , . . . ,x n ) on M, the Poisson tensor P associated with the Poisson 
bracket {•, •} is represented as 

P = Y" piiJL a with pv = {x\ x 1 }. 

/ — ' OXi OXj 

The Jacobi identity is translated into a quadratic differential condition on the matrix 
P 13 known as vanishing Schouten bracket (see, for further details, e.g., [2S1EH]), which 
in local coordinates reads 

n BPfo fjpki BPW 

E pls ^r~ + p -JT- + p -^r- = o, v<> j > (2.i) 

z — ' ox s ox s ox s 

s=l 

We notice that the existence of a Poisson structure does not give any restriction on 
the dimension of M. As an example, the reader should consider the dual q* of a 
finite dimensional Lie algebra g, with Poisson brackets defined, thanks to the natural 
identification q** ~ 0, by 

{f,g}(X) = (X, [df(X),dg(X)\), leg'./.jG C°°(g*). 
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The most familiar instances of Poisson manifolds come from classical mechanics and 
are the cotangent bundles to smooth manifolds. As it is well known (see, e.g., [3]), 
they are equipped a canonical Poisson structure associated with a symplectic, that is, 
closed (indeed exact) and non degenerate, two-form u. A Poisson manifold (M, {-, •}), 
is called symplectic if the Poisson tensor P associated with {-, •} is everywhere invert- 
ible. In such a case, the connection between the Poisson tensor P, the Poisson brackets 
{/, g} and the symplectic two-form u can be described in the following equivalent 
ways: 

a) In local coordinates (x\, . . . , x n ) one has to = J2i j ^ijdxi A dxj, where the matrix 

ujij is the inverse matrix to P % K 

b) It holds the equality 

{f,g} = u>(X f ,X a ) 
where, e.g., Xf = PdF is the Hamiltonian vector field associated with /. 

c) For every vector field Y and every function /, one has: 

Ly(f) = (df 1 Y)=u(Y 1 X f ) 

The local structure of a Poisson manifold can be described as follows (see, e.g., [UHl 
ITT|). Under suitable assumptions on the regularity of the Poisson brackets, one sees 
that an open dense set of M the Poisson tensor has rank r = In. There, M is 
foliated in regular Poisson submanifolds, called generic symplectic leaves, that are 
the (generic) common level sets of k functions Ci, . . . , called Casimir functions of 
P. The distinguished property of the Casimirs of P is that their Poisson bracket with 
any other function on M vanishes, or, equivalently, their differential lies in the kernel 
of P. The dimension of M is related with the integers n and k by dimM = k + In 
As an example, we notice that, in the case of M — sl(n)*, which we can identify 
with sl(n) itself by means, say, of the trace form, the Casimir functions are Cj = 
Tr(X* +1 ), % = 1, . . . , n — 1, and the generic symplectic leaves consist of diagonalizable 
matrices with distinct eigenvalues that can be identified with the direct sum of the 
lower and upper nilpotent subalgebras. 

Let us now come to the definition of bihamiltonian manifold. 

Definition 2.2 A manifold M is called a bihamiltonian manifold if it is endowed 
with two Poisson brackets {f,g}, {f,g}', such that, for any A G M (or A G C if M is 
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complex), the linear combination 

{/, g}' ~ \{f, 9} = (df, (P - \P')dg) (2.2) 

defines a Poisson bracket. This property is known as the compatibility condition 
between the two brackets. 

The expression (12. 2|) will be referred to as pencil of Poisson brackets, and the sum 
Pa = P' — XP pencil of Poisson tensors. The most 'popular' property of bihamiltonian 
manifolds is contained in the following 

Proposition 2.3 Let f and f two functions on a bihamiltonian manifold M , which 
satisfy the characteristic condition Pdf = P'df . Then the Poisson brackets {/, /'} 
and {/, /'}' vanish. 

Proof. It consists of a one-line computation. Let us consider, e. g., {/,/'}: 
{/,/'} = (df,Pdf) = -(df,Pdf) = -(df'P'df) = 0. 

□ 

Definition 2.4 A vector field X that can be written as X = Pdf = P'df is called a 
bihamiltonian vector field. 

Corollary 2.5 Let fk,k G Z be a sequence of (non- constant) functions satisfying 

Pdf l = P'df t+1 . (2.3) 
Then {f u f k } = {ft, f k }' = 0, Vi, k G Z. 

Proof. Using twice equation ()2.3|) and the antisymmetry of the Poisson brackets we 
have, e.g., 

{fi,fk} = (dfi,Pdf k ) = (dfiP'df^) = -(dh-uP'dfi) = 
- (dfk-i, Pdfi+i) = (df i+1 , Ptfk-i) = {fi+iJk-i}- 

Supposing k > i and iterating this procedure k — i times we get {fi,f k } = {fk, fi}- 

□ 



7 



Following ^H], we can state the following 

Amplification 2.6 Let f n ,n = 0, 1, . . . and g n , n = 0, 1, . . . two sequences of (non- 
constant) functions satisfying 

Pdfi = P'df i+l - Pdf = 0, Pd gi = P'dg i+l - Pdg = 0. (2.4) 

Then, along with {f n , f m } = {f n , f m }' = {g n ,g m } = {9n,9m}' = 0, it holds 

{fn, 9m} = {fn, 9m}' = 0, V«, 771 G N. 

The family of vector fields associated with a sequence of functions satisfying the recur- 
sion relations (|2.H|) are customarily said to form a Lenard-Magri sequence. Lenard- 
Magri sequences that start from the null vector field, as in Amplification 12.61 are 
pictorially called anchored Lenard-Magri sequences. Notice that anchored Lenard se- 
quences can occur in bihamiltonian manifold where at east one of the Poisson brackets 
is non-symplectic (indeed, e.g., dfo is a non-trivial element of the kernel of P). We 
can compactly express equations (|2.4jl relative, say, to the sequence /, by considering 
the formal Laurent series /(A) = J2°l fi/X 1 via the equation 

(P' - XP)df(X) = 0. (2.5) 

In analogy with the definition of Casimir of a Poisson bracket, Laurent series satisfying 
(|2.5jl are called Casimirs of the Poisson pencil. The reader should, however, bear 
in mind that while Casimir functions for a single Poisson bracket are, in a sense, 
uninteresting functions, Casimirs of a pencil of Poisson bracket compactly encode 
non-trivial dynamics and constants of the motion. 

If, as it often happens in the applications, inside the family fa we have an element 
/„ satisfying P'df n = 0, we can form a polynomial Casimir of the pencil as 

F(X) = X n f + X n - 1 f 1 + ••• + /„. (2.6) 

Anchored Lenard sequences may give rise to families of integrable systems. Let us see 
how this happens in the case of a In + 1-dimensional manifold endowed with a rank 
2n pencil of Poisson tensors. Let us suppose that we have found a polynomial Casimir 
of the form (|2.6jl . with the n + 1 functions fo, ■ ■ ■ , f n independent. Let S c be a generic 
symplectic leaf of P corresponding to fo = c. The vector fields Xf i} i = 1, . . . , n are 
tangent to S c , are Hamiltonian on S c (w.r.t. P) and the restriction of the functions 
fi, . . . , f n provide n commuting integrals for each of them. In general it holds [TT| I18|: 
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Proposition 2.7 Let (M, P, P') a bihamiltonian manifold of dimension d = 2n + 
k, and let dimiKeriP' — XP)) = k, for generic values of X. Let us suppose that 
H^(X), . . . , H^ k \X) are k polynomial Casimirs of the pencil P\ of the form 

H®(\) = A ni #f + X n ~ l Hf ) +... + H® 

such that the collection of differentials {dHj} 3 ^'"'^ defines a k + n dimensional 
distribution in T*M. Then the vector fields defined by the anchored sequences asso- 
ciated with the 's are integrable (in the ArnoV d-Liouville sense) on the generic 
symplectic leaves of P. 

2.1 Geometry of regular bihamiltonian manifolds and Dar- 
boux - Nijenhuis coordinates 

An important class of bihamiltonian manifold occurs when one element of the Pois- 
son pencil (which without loss of generality we will assume to be P) is everywhere 
invertible, i.e., the Poisson bracket {•, •} associated with P is symplectic. The pos- 
sibility of defining the inverse to one of the Poisson tensors lead us to introduce a 
fundamental object in the bihamiltonian theory of SoV: the Nijenhuis (or Hereditary, 
or Recursion) operator 

N = P'P~\ (2.7) 

(together with its dual N* = P" 1 P'). By definition, N (resp., N*) is an endomor- 
phism of the tangent bundle to M (resp., of the cotangent bundle). As a remarkable 
consequence of the compatibility of P and P' the Nijenhuis torsion of N, defined by 
its action on a pair of vector fields X, Y as 

T(N)(X, Y) = [NX, NY] - N([NX, Y] + [X, NY] - N[X, Y]) (2.8) 

identically vanishes [22] • So, from the classical Frolicher-Nijenhuis theory, we know 
that its eigenspaces are integrable distributions. Such distributions will be the building 
blocks of the bihamiltonian set-up for SoV. 

To explain this, we have to make some remarks and a genericity assumption. It 
can be shown that, owing to the antisymmetry of the Poisson tensors defining N, 
pointwise the eigenspaces of N are even dimensional. Througout the paper, we will 
assume that for generic points m G M, N has the maximal number n = ^dim(M) 
of different eigenvalues Ai, . . . , A n , so that the dimension of the eigenspace relative 
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to any eigenvalue is 2. Otherwise stated, the characteristic polynomial of N is the 
square of the degree n = ~dim(M) minimal polynomial Ajy(A), whose roots are 
pairwise distinct. We will call bihamiltonian manifolds endowed with a Poisson pencil 
with at least one of the elements of the Poisson pencil invertible, and such that the 
eigenvalues of the associated Nijenhuis tensor are maximally distinct, regular biham- 
iltonian manifolds. 

Theorem 2.8 On a regular bihamiltonian manifold there exist a class of coordinates 
(yi,Xi), to be called Darboux-Nijenhuis (DN) coordinates, satisfying the two proper- 
ties: 

(Darboux) They are canonical, that is, {xi,yj} = 5ij, {xi,Xj} = {yi,yj}=0. 

d d 

(Nijenhuis) They diagonalize N , that is, N = J2i ® ^ 7> — ® dxi). 

oyi oxi 

Th proof of this Theorem can be found in |2UEj- Here we will sketch it and discuss 
its meaning. In words, the assertion states that DN coordinates are defined by the 
spectral properties of N, as follows. For all m in the open set U where the eigenvalues 
Aj of iV (which are the same as the eigenvalues of iV*) satisfy Aj 7^ \j, i ^ j, the 
cotangent space T^M admits the decomposition 

T^M = ®™ =1 V mX , dimD m>Ai = 2 (2.9) 

into eigenspaces of A^*. Thanks to the vanishing of the torsion of N, each eigenspace 
Vra^ is locally generated by differentials of pairs of independent functions (fi,gi)- 
This means that the pointwise decomposition (|2.9|) holds (in U' C U) as 

T*M l/v = ®UVx t , 

where V\. is spanned by dfi and dg i: with N*dfc = Xidfi and N*dgi = \dgi. 

Functions whose differential belong to different summands T>\. are in involution 
with respect to the Poisson brackets defined both by P and P' . Indeed, suppose that 
fx and f 2 satisfy N*df x = X x df u N*df 2 = X 2 df 2 , X 1 ^ X 2 . The relation A^* = P^P' 
implies that P'df x = X x Pdf u P'df 2 = X 2 Pdf 2 . So, 



{dh,P'dh) = \ 2 (dh,Pdf 2 ) = A 2 {/ l5 / 2 } 
-{df 2 ,P'dh) = -X x {df 2 ,Pdf x ) = Ai{A,/ 2 } 
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whence the assertion. It is equally straightforward to realize that the only non van- 
ishing Poisson brackets have the form 

{fi,9i} = Fi(f%,9i), {fi,9i}' = F'i{fh9i)i z" = l,...,ra. 

This means that, by quadratures, from the n pairs of functions {/j, gi}%=i n we can 
construct a set of canonical coordinates satisfying the Nijenhuis property of The- 
orem 12.81 Thus the class of coordinates where to frame the bihamiltonian set-up 
for SoV admits a clearcut and simple geometrical description. Admittedly, in the 
general case the computation of DN coordinates requires the integration of the two- 
dimensional distributions T>\. associated with the eigenvalues \ of N*. Fortunately 
enough there are instances (that frequently occur in the applications) in which DN 
coordinates can be found in an easier way. 



2.2 On Darboux-Nijenhuis coordinates 



In this subsection we will briefly discuss conditions and 'recipes' to algebraically find 
and/or characterize sets of Darboux-Nijenhuis coordinates on regular bihamiltonian 
manifolds. A very simplifying instance occurs whenever the eigenvalues \ of N (that 
are, in general, functions of the point m € M) are functionally independent. It holds 
(see, e.g., [23]): 



Proposition 2.9 Let us define h := ^.TrN 1 , i = 1, . . . ,n. In the open set U where 
dli A • • • A • • • dl n 7^ the eigenvalues \ , i = 1, . . . , n are functionally independent, 
XidXi, and so may be used to construct a set of Darboux-Nijenhuis 



satisfy N*dXi 
coordinates. 

Proof. We express the normalized traces Ii of the Nijenhuis tensor N in terms of its 



eigenvalues as kl}~ 


= EL 




Hence dl k = Y%= 




' dh ' 




1 1 




dh 




Ai A2 




dl n 




a" -1 ait 1 ••• 



n-1 





dXi 




dX 2 




dX n 



So we have 



dh A • • ■ A dl n 



X 



Xj)dXi A • • ■ A dX n , 



(2.10) 
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i.e., on the open set where the traces of the powers of the Nijenhuis tensor are func- 
tionally independent, we have that the eigenvalues Aj are different and functionally 
independent. 

To proceed further we need to recall that [23J that the normalized traces Jj of the 
powers of Nijenhuis operator satisfy the recursion relation: 

N*dl k = dl k+1 . (2.11) 

This can be proved as follows. At first one notices that (|2.11j) is equivalent to the 
relation 

L N x{h) = L x (I k+1 ), V vector field X, 

as it can be easily seen evaluating the equality of one-forms on a generic vector field 
X. Thanks to the Leibnitz property of the Lie derivative and the ciclicity of the trace, 
we see that 

L NX (I k ) = Tr(L NX (N) ■ N"- 1 ) and L x (I k+1 ) = Tr(L x (N) ■ N k ). (2.12) 

Since the vanishing of the Nijenhuis torsion of iV implies that Ln X (N) = N ■ 
L X (N), WX the validity of (EHTJ) is proved.. 

We now express the relations ()2.11|) in terms of the eigenvalues Aj 



1 1 


1 




' N*dX x 


— \1dX1 




' " 


Ai A2 


X n 




N*d\ 2 


- A 2 dA 2 







\n— 1 \n— 1 
A l A 2 


\n— 1 




N*d\ n 


— X n d\ n 








Since the Vandermond matrix in the LHS of this equation is, by assumption, 
invertible, we conclude that, for i — 1, . . . , n, it holds N*dXi = XidXi. 

□ 

This proposition can be rephrazed saying that "half of" the DN coordinates are alge- 
braically provided by the Nijenhuis tensor itself. The remaining "half" yi, . . . , y n can 
always be found by quadratures. Actually, there is a condition leading to the alge- 
braic solution of this problem too. To elucidate this, the following two considerations 
are crucial. 
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First of all one remarks that, thanks to the fact that, for DN coordinates, if i ^ j 
one has {Xi,Vj} = one can replace the n(n — l)/2 equations {\i,yj} = Sij with the 
n equations 

{Ai H h X n ,yj} = 1, j = 1, • 

that do not involve the individual coordinates Aj but only their sum J^" =1 A, = 7i, 
and the Hamiltonian vector field 

The second argument goes as follows. Let us consider the distinguished functions Ik 
introduced in propositior J2.9l and trade them for the coefficients pi of the minimal 
polynomial 

A N (X) = X n -p 1 X n - 1 -p 2 X n ~ 2 Pn 

of N. The functions pk and Ik are related by triangular Newton formulas 

1 2 1 , 
h=Pi, h=P2 + jzPi, h=ps+ P2P1 + qPi; 

2 ! ! i 3 (2-14) 
7 4 = Pi + pm + vlvi + + J 5 = P5 + ■ ■ ■ ■ 

As a consequence of the recursion relations (J2.11j) . it can be easily shown that the 
pi's satisfy the 'Frobenius' recursion relations 

N*dpi = dpi+i + pidp!, with p n+1 = 0. (2.15) 

We can compactly write these relations as a single relation for the polynomial Ajy(A); 
indeed, a straightforward computation shows that they are equivalent to 

N*dA N (X) = XdA N (X) + A N (X)d Pl . (2.16) 

Actually relations of this kind are very important for our purposes. Indeed, in ^H] 
we proved the following proposition: 

Proposition 2.10 Let $(A) a smooth function defined on the manifold M , depending 
on an additional parameter X. Suppose that there exists a one-form a$ such that 

N*d$(X) = Arf$(A) + A N (X)a<s, . (2.17) 

Then, the n functions $j obtained evaluating the "generating" function $(A) for X = 
Aj, i = 1, . . . , n are Nijenhuis functions, that is, they satisfy N*d$i = Ajd$j. 
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Definition 2.11 We will call a generating function $(A) satisfying equation \2.1l ) 
a Nijenhuis functions generator. 

The relevance of the notion of Nijenhuis functions generator in the search for DN 
coordinates stems from the fact that Nijenhuis functions generators form an algebra 
jV(M), which is closed under the action of the vector field Y = —Pdl\. In this way, 
knowing a set of Nijenhuis functions generators, we can obtain further elements of 
the algebra M{M) by repeated applications of the vector field Y. Clearly, in such an 
extended algebra, the characteristic equation 

M*(A)) = 1 + A JV (A)a v , 

may be easier to be solved, thus yielding the missing Darboux-Nijenhuis coordinates 
coordinates yi as yi = \I/(Aj). We will see an instance of this situation in the Example 
of Section HJ 

For an analysis of Darboux-Nijenhuis coordinates within the theory of multi- 
hamilronian structure on loop algebras, see [HI El HHj 



3 The Separability Conditions 

As we have briefly recalled in Section on a bihamiltonian manifold one is usually 
lead to consider bihamiltonian vector fields, that is vector fields X admitting the 
twofold Hamiltonian representation X = Pdf = P'dg. Let us now suppose that 
(M, P, P') be a regular bihamiltonian manifold of dimension 2n, and that we were 
able to construct, by means of the Lenard-Magri iteration procedure a sequence of 
functions Hi, H 2 , . . . satisfying P'dHi = PdH i+ i. Let us also suppose that the first n 
of them be functionally independent. Then one easily shows that the all the further 
Hamiltonians H n+ i, . . . are functionally dependent from the first n. Indeed this follows 
from the fact that a regular Poisson manifold of dimension In cannot have more than 
n mutually commuting independent functions. 

This means that, if we consider the Hamiltonian H n+ i, there must be a relation 
of the form 

^(Hi,...,H n ,H n+ i)=0, with ^, Hn+1 = + 0. (3.1) 

relating it with the independent Hamiltonians Hi, % = 1, . . . , n. 
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Actually, the case of Hi = ij = ^.TrN 1 is an instance of this situation. In fact, 
since by the Cayley-Hamilton theorem N annihilates its minimal polynomial, we have 
N n — YH=iViN n ~ l = 0, yielding the relation 

n 

I n+ x/2(n + 1) - J2 2 ( n ~ 1 + iWn-i+i = °- 
i=i 

Differentiating equation (jHHj) we see that, along with P'dHi = PdH i+ i, i — 1, . . . , n—1 
it holds: 

P'dH n = PdH n+1 = ( V ^-PdHi) , (3.2) 

that is, the vector field X n+ i = PdH n+ i = P'dH n is a linear combination of the 
vector fields X 1 = PdH 1} ...X n = PdH n . 

This innocent looking observation is the clue for the bihamiltonian theory of SoV. 
Let {Hi, H 2 , ■ ■ ■ H n } any integrable system on M, that is that is, the if^s are mutually 
commuting (w.r.t.P) independent functions. We can construct a n-dimensional dis- 
tribution, namely the distribution T>n spanned by the n mutually commuting vector 
fields Xi = PdHi. This is nothing but the very classical tangent distribution to the 
invariant tori of the Liouville Arnol'd theory of integrable systems. Since M comes 
equipped with a second Poisson tensor P' we can as well consider the distribution 
T>' H generated by the Hamiltonians Hi under the action of P' , that is, generated by 
the vector fields X[ = P'dHi. It holds: 

Theorem 3.1 Let {Hi, . . . , H' n } define, as explained above, an integrable system on a 
regular bihamiltonian manifold (M, P, P'). The Hamilton- Jacobi equations associated 
with any of the Hamiltonians Hi are separable in the DN coordinates x%)i defined by 
N = P'P' 1 if and only if the distribution T>' H is contained in T>h, or, equivalently, if 
the distribution T>h is invariant along N . 

Proof. We will first prove the equivalence of the invariance of T>h under N and the 
inclusion T>' H C T>h- To say that T>' H is contained in T>h is tantamount to saying that 
there exists a matrix Fij, whose entries are, in general, functions defined on M such 
that, for i — 1, . . . , n it holds: 

X[ = P'dH = J2 F ij pdH i = Yl Fi i X r ( 3 - 3 ) 

3 j 
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Writing P' = NP we can translate these equalities into NXi = FijXj for all 
i = 1, . . . , n. 

The full proof of the fact that the invariance of T>h insures separability in DN 
coordinates can be found in ^Hj- It goes as follows. 

At first we notice that the translation in terms of the codistribution D* H generated 
by the differentials of the Hamiltonians Hi of the invariance condition for T>h is the 
invariance condition N*T>* H C T>* H . This can be easily seen applying to ()3.3j) the 
operator P" 1 , to get N*dHi = Y,j F ij dH j- 

Since all the Poisson brackets {Hi,Hj} vanish and M is a regular bihamiltonian 
manifold, the matrix F defined by ()3.3)1 can be show to have simple eigenvalues, that 
coincide with eigenvalues A, of N. So there exists a matrix S satisfying 

SF = AS, where A = diag(Ai, . . . , A n ). 

If we define the n one forms 6i = ^\ SijdHj, we get: 

N*9i = SylTdHj = SijF^dH k = ^ AA ; .V,,r///, = (3.4) 

j 3,k j,k 

This means that 9i is an eigenvector of A^* relative to A^; hence there must be functions 
Fi, Gi such that 

2J SijdHj = Fidxi + Gidyi (3.5) 
j 

whence the existence of a separation relation yi] Hi, . . . , Hh) for all i — 1, . . . , n. 

The converse statement can be trivially proved. 

□ 

We would like to stress thatthe separability condition of Theorem 13.11 is a tensorial 
one. That is, given a regular bihamiltonian manifold (M, P, P') this separability 
criterion can be checked in any system of coordinates, without the a-priori calculation 
of the DN coordinates themselves. Notice, also, that the validity of the statement 
does not (as it should be!) depend on the choice of mutually commuting integrals 
{Hi, . . . ,H n }. That is, if we consider a "change of coordinates in the space of the 
actions", that is we trade the if;'s for another complete set of integrals of the motion 
Ki = Ki(Hi, . . . , H n ), then the separability of the new Hamiltonians Ki will hold if 
and only if the separability of the original ones holds. Indeed, the dual distributions 
generated by the H^s and the K^s coincide. 
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A second remark is important and deserves to be explicitly spelled out. Although 
we have started our discussion considering the case of a family of bihamiltonian 
vector fields, that is the case of Lenard-Magri sequences, the hypotheses of Theorem 
13.11 concern only the relations of the distributions generated respectively under the 
action of P and P' by the Hamiltonians Hi, without any mention of the fact that 
the generators of the distribution be bihamiltonian vector fields. Thus, although it 
might seem a somewhat odd statement, the vector fields that are separable by means 
of the bihamiltonian approach are not necessarily bihamiltonian vector fields!. It is 
also important to notice that it is not only a matter of choice of generators. Indeed, in 
|2Zj it has been shown that the only bihamiltonian vector fields on a regular bihamil- 
tonian manifold turn out to be associated with separated functions of the eigenvalues 
of N, that is functions of the form H = Y^i=i fi{^i)- This means that, in such a case, 
the distribution T>h coincides with that generated by the distinguished functions U. 
However, this is by far a very special example, that is, the range of applicability of 
the method is much wider than that, as it as already been quite widely shown in the 
literature. 

The separation condition of Theorem 13 .11 is based on the analysis of the behaviour 
of the characteristic distribution associated with an integrable system under the Nijen- 
huis tensor N . An equivalent criterion, based on the analysis of the Poisson brackets 
associated with the tensor P' can be formulated as follows. 

Theorem 3.2 Let {H = Hi, H2, . . . H n } an integrable system defined on regular 
bihamiltonian manifold (M, P, P') . The Hamiltonians Hi are separable in the DN 
coordinates defined by N = P'P^ 1 if and only if, along with the commutation rela- 
tions {Hi,Hj} = there hold also 

{Hi, Hj}' = (dHi, P'dHj) = 0, for i,j = l,...,n. (3.6) 

Proof. The key formula is the relation between P,P' and N*. Indeed, suppose that 
T>* H be invariant along N*. Then: 

{Hi,HjY = (dHi, P'dHj) = (dHi,NPdHj) 

(WdH^PdHj) = Y, F ik(dH k ,PdH 3 ) = Y,F lk {H k ,Hi] = 0, 

k k 

which, in view of Theorem 13.11 proves the statement in one direction. Now, let us 
suppose that (|3.6|) holds. Then, for every ij — 1, . . . , n we have: 

= {Hi,Hj}' = (dH h P'dHj) = (dH h NPdHj) = (N*dH h PdHj), 
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meaning that, for all i = 1, . . . , n, the one-form N*dHi belongs to the annihilator 
(w.r.t. P) of the distribution T>h- Since such an annihilator coincides with T** H , this 
means that, for alH = 1, . . . , n, N*dHi 6 T>* H . 

□ 

This results leads to the following, somewhat daring, comparison. The Liouville- 
Arnol'd theorem on finite dimensional integrable Hamiltonian systems says that the 
geometrical structure underlying integrability of a Hamiltonian vector field defined 
on a symplectic manifold M,us is a Lagrangian foliation of M. We can rephrase 
the content of Theorem 13.21 saying that the geometrical structure underlying the 
separability of a system defined on a regular bihamiltonian manifold(M, P, P') is a 
bilagrangian foliation of M. 

We end our presentation of the bihamiltonian set-up for SoV with the following 
remark. Theorem 13.11 concerns only the existence of the separation relations. In 
principle, one could try to find these relations in concrete examples by actually di- 
agonalizing the matrix F, and explicitly finding and integrating the relations ()3.5j) . 
However, there is a very simple tensorial criterion which can be used to determine 
the functional form of the separation relations yi, Hi, . . . , H h ), whose proof can 

be found in [To] . 

Proposition 3.3 Let {Hi, . . . , H n } be an integrable system defined on a regular biham- 
iltonian manifold, which is separable in the Darboux-Nijenhuis coordinates associated 
with N = P _1 P'. Consider the matrix P^ fulfilling the relations hS.^l . Then the 
separation relations are be affine in the Hamiltonians H , that is of the form 

$i{xi, yi, Hi,..., H h ) = 22 s ij{ x h Vi)Hj + Ui(xi, yi), (3.7) 

j 

if and only if the matrix F satisfies the relation N*dFij = F^dF^. 

The matrix S on ()3.7j) can be shown to be a suitably normalized matrix of eigenvectors 
of the matrix P. Its characteristic property is that, as expressed in the equation, the 
entries Sij of the i-ih row depend only on the pair (x^yi) of Darboux-Nijenhuis 
coordinates. For this reason it can be called a Stackel matrix. 
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4 Example: a generalized Toda Lattice 



In this final Section we will apply the general scheme outlined in the previous Sections 
to a specific model, with the aim of showing how the recipes discussed so far from 
a theoretical standpoint can be concretely applied. We will study a generalization 
of the four site Toda lattice, to be termed Toda\ model. This system is a member 
of a family introduced in [21] as reductions of the discrete KP hierarchy. It can be 
described as follows. We consider on M = C 12 , endowed with global coordinates 
{bi, di, Ci}j=i..4 the Hamiltonian 



h gt = + bl + h + b\) - (ai + a 2 + a 3 + a 4 ), 
and the linear Poisson tensor given by the matrix 



(4.1) 



P 











— CLi 





a 4 


Ax 


Cx " 




a% 


— a 2 





c 2 





with A x = 



















a 2 


-a 3 












o 3 — 04 



(4.2) 





" -c x 





c 3 


" 







-c 2 





c 4 


Cx = 












Cx 





-c 3 










c 2 





-c 4 . 



c 5 



' 


~Cx 





c 4 


Cl 





-c 2 








c 2 





-c 3 


. ~ca 





c 3 






where we denoted by the 4x4 matrix with vanishing entries.. 
Using (here and in the sequel) the cyclic identifications aj +4 = aj, b i+i = 
the Hamiltonian vector field Xh gt = PdHcT can be written as follows 



bi, Ci+4 



h 




a i-l — a i 


hi 




ai{b i+ x ~ 


Ci _ 




_ Q(6i_2 - 



h) + c,i 



,4, 



(4.3) 



The expert reader will notice that Hgt coincides with the Hamiltonian of the periodic 
four-site Toda lattice, written in Flaschka coordinates bi = Pi,a>i = exp(gj — q- t+ x)- 
Indeed, on the hyperplane Mt — C 8 defined by q = 0, i = 1, . . . 4 the vector field 
Xh g t defines the periodic Toda flow. 



19 



Proposition 4.1 The Hamiltonian vector field Xh gt admits the Lax representation 
L(n) = [L(/i),$], where 



L{n) 







Ol 

M 


b 2 


ci 

w 


22. 
1' 




C2 



C2 Q3 h 



£4 

,/2 









C3 
/? 


Q 4 
/' 


£1 

/' 








C4 


ci 
J? 


A' 











C2 
/? 


03. 
l" 






(4.4) 



The bihamiltonian aspects of this system have been discussed in (2H1- I n particular, 
it has been noticed that on M there exists a second Hamiltonian structure for the 
vector field Xh gt - Namely one considers the bivector P' having the following form: 



P' 



A 2 


Bx 


c 3 


-Br 


A 2 


c 4 




-CJ 


A 4 



where C3 



" -61C1 













-b 2 c 2 





6 2 c 4 







-&3C3 








C264 





— 64C4 



(4.5) 





" 







— a 4 




-a 1 










A 2 = 















— a 2 









a 4 





-a 3 






Ac 





&2C1 + a x a 2 


— &1C4 — O1O4 
-aiCi — aiC2 CI1C3 
cia 2 — a 2 c 2 — a 2 c 3 
cia 3 c 2 a 3 -a 3 C3 
— cia 4 C2a 4 C3CI4 



Bx 



-6 2 ci — aia 2 


63 c 2 + a 2 a 3 


aic 4 

a 2 c 4 
-a 3 c 4 
— a 4 c 4 



-6iai Ci — C3 6ia 4 
6 2 ai -&2«2 c 2 -c 4 
-ci 6 3 a 2 -&3«3 c 3 
c 4 — C2 6 4 a3 — & 4 a 4 
&ic 4 + aia 4 

-63C2 - a 2 a 3 

— 6 4 C3 — 0304 

&4C3 + a 3 a 4 

— C1C2 C1C4 
cic 2 -C2C3 
c 2 c 3 -C3C4 

— C1C4 C 3 C4 



20 



It can be easily noticed that Xh gt = P'd(— £) i=1 h)- More in general, we have the 
following 

Proposition 4.2 The pencil P' — XP is a pencil of Poisson brackets. The rank of 
the generic element of the pencil is eight. The characteristic polynomial R(X, p) = 
DetiXl — L(/i)) ; written in terms of p = p 4 can be expanded as: 

R(X, p) = A 4 - p + H(X) + (K(X) - X 2 J 1 )/p + J 2 /p 2 . (4.6) 

The functions J% and J 2 are common Casimirs of P and P' . The polynomials 
H(X), K(X) are polynomial Casimirs of the pencil P\ — P' — XP. They have the 
form 

H(X) = X 3 H - X 2 Hi + XH 2 - H 3 , K(X) = K X + K x . (4.7) 

Explicitly, J\ = Cic 3 + c 2 c 4 , J 2 = C1C2C3C4, while the coefficients of H(X) and K(X) are 
given by: 

4 4 4 4 

H = ^2 b i, Hi = ^2 hbj + ^ a " H z = ^2 ( Q + hi ^ ai + 1 + fll + 2 ) + bib i+ ib i+2 ) 

i=l i>j=l i=l i=l 

4 

H 3 = hc i+ i + aia 3 + a 2 a 4 + cubic and quartic terms; 
i=i 

4 4 

K = bjCj-iCj+t - Cia-i-ia i+2l K\ = OiC 1+1 c 1+2 + quartic terms 
i=i i=i 

One can show that the eight functions Ho, Hi, H 2 , H3, Kq, Ki, Ji, J 2 are functionally 
independent and, thanks to the fact that they fill in Lenard sequences, are mutually in 
involution. The kernel of P is generated (at generic points m G M) by the differentials 
of the four functions Ho, Kq, Ji, J 2 . Hence, on the 8 dimensional manifold S K defined 
by the equations H = Ki, K = K 2 , Ji = ^3, J2 = that is, the generic symplectic 
leaf of P, the vector field Xh gt is completely integrable. To realize this we simply 
have to notice that H GT can be expressed as \Hq — Hi, and apply the properties of 
anchored Lenard Magri sequences collected in Proposition 12.71 

4.1 Separation of Variables 

We will now show how to apply the ideas and recipes if the bihamiltonian set-up for 
SoV to the Toda 4 model introduced above. The first problem to deal with is that the 
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Poisson tensor P' does not restrict to S K , but must be projected. This can be done 
as follows (see [311311111211 for details and the geometric background), by means of 
a kind of Dirac reduction process. 

d d 

We consider the vector fields Z\ = —7-— and Z 2 = — — , notice that the matrix 

004 ' 

L Zl (H ) L Zl {K x ) 
L Z2 (H ) L Z ,{K X \ 



G = 

is invertible, and form the bivector 




2 

A = ^2 [G^jZi A X{, where X\ = P'dH , X\ = P'dK . (4i 



Lemma 1 The modified bivector Q = P' — A defines a Poisson bracket, compatible 
with P; Q restricts to S K . 

Proof. The proof of the fact that Q\ = Q — XP is a Poisson pencil follows (see, e.g., 
[To]), from the equalities: 

L Zl P = 0, L Zl P' = W\ NZ X - c 3 ^- A Z 2 

3 (4 9) 

L Z2 P =(—-—) A Z 2 , L Z2 P' = (6 4 — + — ) A Z x + Wi A Z 2 . 
obi oai ' aa 4 obi 

d d d 

with W\ = (a 3 - a 4 - hc 2 - c 4 — ) and 

oas 004 oc 2 0C4 

ttt2 n ® 7 \ 

W 2 = ( b 4-^T - h l^T ~ a 3a C l^— + C 2TT~ + C 3 ~ C 4 — ), 

ao 4 obi oa 3 oa^ Oci oc 2 oc 3 oc^ 
as well as from the fact that 

QdH = QdK = dJi = QdJ 2 = 0. (4.10) 

To show that (|4.9|) holds true is a matter of an explicit computation, while (j4.10|) 
follows from the definition of Q. In fact, the last two equations hold since J\ and J 2 
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are Casimirs of P' invariant under Z\ and Z 2 . For, e.g., H one computes 

2 

QdH = P'dH - AdH = X{ - {[Gij]L Zj (H )) ■ X\ 

2 

= x l - E ([^1^0 ■ x i = *J - E ki ■ x i = °. 

ij=l i 

where the second equality follows from the fact that all the functions Hi, K a , J a are 
in involution w.r.t. P. 

□ 

Thanks to the above lemma, the generic symplectic leaf iS re is endowed with the struc- 
ture of a regular bihamiltonian manifold. It is easy to show that the non trivial 
Hamiltonians Hi, H 2 , H 3 , Ki (more precisely, the restriction to iS^of these Hamiltoni- 
ans) satisfy the hypothesis of Theorem 13.11 w.r.t the (restriction to S K ) of the pencil 
Q — XP. Indeed we have: 

2 2 

QdH, = P'dH % - [G-}]Zi A X{{dHi) = PdH i+1 - £ [G^L^H^Xi 

i,j=l i,j=l 

(where we understand H± = 0) and 

2 2 

QdK 2 = P'dKi - J2 [G7j]Zi A X{{dK 2 ) = Y [G^)L Zi {dK x )X{ 

i,j=l i,j=l 

So we proved that, for generic values of the Casimirs i = 1, . . . , 4, the system 
obtained by restriction of the Todag flows on S^is separable in the DN coordinates 
associated with the restriction to S^of pencil Q — XP. To finish our job we finally 
have to: 

a) explicitly define the DN coordinates; 

b) find the separation relations. 

To solve the first problem, we will use the tools briefly described in Section 12.21 
We rely on a result of [TE], as well as on explicit computations, to state the following 
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Proposition 4.3 Let Q(X) be the matrix 



L Zl H(X) L Z2 H{X) 
L Zl K(X) L Z2 K(X) 



(4.11) 



The roots of the degree 4 polynomial Det{Q) are the roots of the minimal polynomial 
A(A) = A 4 — Y^^iPiX 4 ' 1 of the Nijenhuis tensor N = P~ l Q associated with the 
regular Poisson pencil Q\. The coefficients p% are well defined functions on the generic 
symplectic leaf S K , and are functionally independent. Furthermore, the ratios p(X) = 
—Qi,%l$Xfr c(X) = —Qi,xjQx,\ are Nijenhuis function generators. 

Thus, one half of the Darboux-Nijenhuis coordinates will be given by the roots 
of Det(£?). To find the remaining half we consider vector field Y = —Pdp\ whose 
role has been discussed in Section 12.21 Since an explicit computations shows that 
L Y log(p(A)) = 1, hence we can state the following 

Proposition 4.4 A set of Darboux-Nijenhuis coordinates for the restriction to the 
generic symplectic leaf S K of the Toda\ flows are given by the four roots Aj of Det(Q)(X) 
and by the values \ii of the function log(p(A)) for X = Aj, where 



P(A) 



— cia 3 - aic 2 ) A + c 2 ai& 3 - aia 2 a 3 + cib 2 a 3 + c x c 2 
CiC 3 A + aia 2 c 3 - ci& 2 c 3 



To find the separation relations we reconsider the Lax matrix f!4.4j) . and notice that 
the pairs, e.g., Aj,p(Aj) are common solutions to the system 



PQ\,\ + 021 
P01,2 + Q22 




0. 



(4.12) 



Since the Lie derivative of the matrix C(X, //) = Al — L(n) along the vector fields Z { 
are given by: 
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L*(£(A,/i)) = 














, L Zl (£(A,/i)) = 

























1 

















,-1 -1 
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we see that the solutions of the system (|4.12|) are nothing but the solutions of the 
equations 

[£(A,/i) v ] 44 = 0, [£(A,/i) v ] 4 i = 0, 

C(X, /i) v being the classical adjoint to £(A, fi) = Xl—L(/j,). Using standard arguments 
of linear algebra, and the results collected in Proposition 14 A\ we can state 

Proposition 4.5 The separation relations connecting pairs of Darboux-Nijenhuis co- 
ordinates Xi, Hi, the Hamiltonians Hi, H 2 , H%, K 2 and the Casimirs H$, K\, K 2 , Jq are, 
on the generic symplectic leaf S K , given by the evaluation of the characteristic poly- 
nomial Det(C(X, //)) in X = Xi, h = Hi = log(p(A,)). 

We notice that, a posteriori, the separation coordinates for the Toda| system fall in 
the class described in, e.g., (331 EJ EH EIJ 13 123 • Namely, the DN coordinates that 
separate the Todag system are Algebro-geometrical Darboux coordinates associated 
with the spectral curve ()4.6|) . and fulfill the so-called Sklyanin's 'magical recipe'. 

As a final remark, in connection with the discussion on the relation between 
the bihamiltonian property of an integrable vector field and the separability of the 
associated HJ equations of Sectional we notice that the Hamiltonians Hi, H 2 , H%, K 2 
are functionally independent from the coefficients of the minimal polynomial of the 
Nijenhuis tensor obtained from Q — XP. So, this is a further instance of a system 
which is not bihamiltonian on a regular manifold, but turns out to be separable via 
the bihamiltonian method of SoV. 

4.2 A remarkable subsystem: the open Toda| system 

In this last subsection we will discuss a remarkable reduction of the periodic Toda| 
system, leading to the corresponding generalization of the open (or non-periodic) one. 
In the manifold M ~ C 12 we consider the nine-dimensional submanifold M defined 
by the equations: 

a 4 = c 3 = c 4 = 0. (4-13) 

One can easily verify that the restriction A^ gt to Mq of the vector field Xh gt is 
tangent to M . Also, the tensor P can be restricted to M ; indeed, the expression of 
its restriction P w.r.t. the natural coordinates {bi, . . . , 64, a%, . . . , a 3 , c 4 , c 2 } of M is 
obtained from ()4.2|) simply by removing the 9th, 11th, 12th rows and columns. 
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One can easily check that X^ gt 



H, 



GT,0 



5« 



P^dHcTfi with 

- (d + a 2 + o 3 ), 



6^ + 6 3 + 6; 



(4.14) 



and recognize that this function is the Hamiltonian of the open Toda lattice. Also, a 
Lax pair for X^j gt is 
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ai 

£L 

n2 



(±2 

C2 
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(4.15) 



It should be clear from the form of the Lax pair that the vector field X^ gt on 
Mo is an extension of the standard open Toda lattice towards the so-called full open 
Toda lattice, which is a system describing a flow on the lower Borel subgroup of 
sl(N). The integrability of the full open Toda lattice was established in jH] (see, also, 
[12 ). The idea was to complement the integrals of the motion coming from the Lax 
representation with additional integrals obtained by means of the 'chopping method', 
whithin the group-theoretical point of view. 

The need to supply the standard results of the Lax theory with further methods 
should be clear from the following considerations. The only Casimir function of Pq 
is h = Y^i=i b%- Hence, its symplectic leaves C (M ) are the eight dimensional 
manifolds defined by ho = £, and X^ gt can be seen as a Hamiltonian system with 
four degrees of freedom. 

The characteristic polynomial of the matrix L is 



Det(Al - L (aO) = -/i 4 + A 4 - h X 3 + M 2 - M + h 3 , 
that is, it provides us with only three non trivial Hamiltonians 



(4.16) 



4 

£ 

i>j=i 



hbj + y^Qj, h 2 = ^2 b i b jh + a *( &i + 2 + b i+s) + ci + c 2 



i=l i>j>k=l i=l 

h-i = bia 2 b 4 + ai6 3 6 4 + 6162CI3 + 6ic 2 + 0103 + C164 + 6i6 2 6 3 6 4 . 



(4.17) 



We will now show how the the tools we previously introduced can be used to 
geometrically prove the complete integrability of such a system and, moreover, yield 
the existence of an additional integral of the motion . 
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The main property is that, along with P, the tensor Q resticts to M$. This can 
be proven as follows: one checks by direct inspection that this is true for P'\ the 
assertion follows from the fact that the vector field X\ (defined in (|4.8jl ) vanishes on 
Mo, while Z\ and X\, which coincides with Xqt, are tangent to Mq at the points of 
M . 

Furthermore, we add two observations. The first one concerns the restriction Q 
to M of the matrix Q. It has the form 



Go 



^11 




A 2 - (62 + 63) A + b 2 b 3 + a 2 
Jcia 3 + aic 2 ) A - aia 2 a 3 + c x c 2 + Ci& 2 a 3 + aic 2 6 3 



(4.18) 



with C/jj = A 3 — 71"! A 2 — 7r 2 A— 7r 3 a degree three polynomial. Hence its determinant (that 
is, the minimal polynomial of the Nijenhuis tensor N induced by the pair Q Q — XP 
on S%) factors as Q\\Q 22 - 

The second observation consists in the fact that the three surviving Hamiltonians 
hi,h 2 , h$ given by (j4.17|) satisfy the conditions: 



Q dh i = Y, F ? j p o<thi, with** 



TTt 1 
7T 2 1 
7T 3 



(4.19) 



We notice that the functions 7T,-,i = 1, 



3 and the root 



A, 



-aia 2 a 3 + c x c 2 + aic 2 6 3 + Ci& 2 a 3 



(4.20) 



Cia 3 + a x c 2 

of (?2 2 are s ^ m functionally independent and hence (generically) different on S^. 

Lemma 2 Let a any function satisfying Q$da = X^P^da. Under the above hypotheses, 
the brackets {a,hi}p and {a ) h i \q vanish. 

Proof. Evaluating both sides of Qoda = X^P^da on the differentials (dhi,dh 2 ,dh 3 ), 
and switching the action of the Poisson tensors on the dh^s, we get 

(da, Q dhi) = X±(da, P dhi),i = 1, ... ,3 

Inserting (|4.19|) we get the equation Yu^=i {Ffj — A 4 5jj) (da, P dhj) = 0. Since A 4 is 
not an eigenvalue of F^, the Lemma is proved. 
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□ 



So a fourth integral of the motion, that commutes with the Hamiltonian Hq T for the 
open Todaf lattice is given indeed by the distinguished root A 4 of equation ()4.20|) ; 
this constructively proves the integrability of the system. 

Finally, we notice that this method proves the existence of a fifth integral of the 
motion. Indeed, we know that, along with A4, there must exist another independent 
function /Z4, satisfying the hypotheses of Lemma El and functionally independent of 
A4 and of the h^s. In such a comparatively low dimensional case, such a function can 
be explicitly found to be 

c 2 (ai& 2 ci - a 2 a 1 2 - c x 2 - Ci& 3 ai) 
ci (cia 3 + a x c 2 ) (A 4 - 7TiA 4 - 7r 2 A 4 - 7r 3 ) 
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